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We construct a random matrix model for the bijection ^ be- 
tween clas- sical and free infinitely divisible distributions: for every 
d> 1, we associate in a quite natural way to each *-infinitely divisible 
distribution n a distribution ¥'^ on the space of d x d Hermitian ma- 
P^ , trices such that P^ * P^ = PJ^*". The spectral distribution of a random 

Q, ■ matrix with distribution P^' converges in probability to *(/i) when 

d tends to +00. It gives, among other things, a new proof of the 
almost sure convergence of the spectral distribution of a matrix of 
C^ ' the GUE and a projection model for the Marchenko-Pastur distri- 

bution. In an analogous way, for every d > 1, we associate to each 
^-infinitely divisible distribution n, a distribution L^ on the space of 
complex (non-Hermitian) d x d random matrices. If /j is symmetric, 
Cn ' the symmetrization of the spectral distribution of |Mtj|, when Ma is 

►^ , L^-distributed, converges in probability to ^(/i)- 

00 

^D ' Introduction. Free convolution ffl, defined in Bercovici and Voiculescu 

(1993), is a binary operation on tiie set of probability measures on the real 
^^ I line, arising from free probability theory (/x ffl z^ is the distribution oi X + Y 

(^ ■ when X,Y are free and have distributions /i, z^). It is associative, commu- 

tative and continuous with respect to the weak convergence. A probability 
measure // on M is said to be ffl-infinitely divisible if for every n > 1, there 
exists a probability measure /u„ on M such that //^" equals to /_i. 

It is shown in Bercovici, Pata and Biane (1999) that there exists an home- 
omorphism ^ from the set of *-infinitely divisible distributions to the set of 
ffl-infinitely divisible distributions which associates to every classical (resp. 
^ • free) limit theorem a free (resp. classical) analogue. Indeed, for every *- 

infinitely divisible distribution //, for every sequence (//„) of probability mea- 
sures, for every sequence (kn) of integers tending to infinity, the sequence 
/i*'^" tends to fi if and only if the sequence /U^'^" tends to ^(^). 



Received February 2004; revised April 2004. 

AMS 2000 subject classifications. Primary 15A52, 46L54; secondary 60E07, 60F05. 
Key words and phrases. Random matrices, free probability, asymptotic freeness, free 
convolution, Marchenko-Pastur distribution, infinitely divisible distributions. 



This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Probability, 
2005, Vol. 33, No. 3, 1134-1170. This reprint differs from the original in 
pagination and typographic detail. 

i 



2 F. BENAYCH-GEORGES 

The proofs in Bercovici, Pata and Biane (1999) rely on integral transfor- 
mations and complex analysis. We will, in this article, construct a niatricial 
model for the ffl-infinitely divisible distributions, and present in a more pal- 
pable way the bijection ^. 

Let fj. be an *-infinitely divisible distribution. Let (^.„) be a sequence 
of probability measures and (kn) a sequence of integers which tends to in- 
finity such that the sequence fi"^" tends weakly to //. Let, ioi d> 1 and 
n > 1, Q|^" (resp. K^") be the distribution of f7diag(X„,i, . . .,Xn,d)U* [resp. 
of U diag{Xn,i, ■ ■ ■ , Xn,d)V], where U,V are independent unitary Haar dis- 
tributed random matrices independent of the i.i.d. random variables Xn^i, • • • , Xn^d 
with distribution /i„. We will prove, in Section 3 (resp. Section 7.1), that 
the sequence ((Q(J")*'^") [resp. ((K^")*'^")] converges weakly to a probabil- 
ity measure P^ (resp. L^). The main results of this article are the following 
ones: the spectral distribution of a random matrix with distribution P^ con- 
verges in probability to ^(^) when d tends to infinity, and so does the 
symmetrization of the spectral distribution of \Md\ when Md is distributed 
according to L^. So we have constructed matrix models which go from *- 
infinitely divisible distributions to ffl-infinitely divisible distributions when 
the dimension goes from one to infinity. What is more, for all *-infinitely 
divisible distributions fi, v and ah d, P^'' = P^ * P^ and L^'' = L^ * L^ . This 
property (and the fact that all formulas depend analytically on d, so could 
be extended to noninteger d) opens the perspective of a continuum between 
the classical convolution * and the free convolution ffl for infinitely divisi- 
ble mesures [M. Anshelevich has already constructed such a continuum in 
Anshelevich (2001), but the model we present here does not interpolate his 
construction]. T. Cabanal-Duvillard, in Cabanal-Duvillard (2004), has stud- 
ied at the same time as the author the distributions P^, and has proved the 
same result, but with different methods (processes, measure concentration, 
integral transforms). 

At last, in the case where [i is the standard normal distribution, ^(/u) 
is the semi-circle distribution with center zero and radius two, and the dis- 
tribution P^ is closely related to the one of the GUE, so that the conver- 
gence of the spectral distribution of a matrix with distribution P^ implies 
Wigner's result. Likewise, the distribution L^ is the one of a matrix with 
independent Gaussian entries, and we have a new proof of the convergence 
of the spectral distribution of the Wishart matrix with parameter 1 to the 
Marchenko-Pastur distribution. 

In the same way, in the case where \i is the classical Poisson distribution, 
this result allows us to see the Marchenko-Pastur distribution as the limit 
spectral distribution of a sum of independent rank-one projections. 

The text is organized as follows. In Section 1 we recall a few results about 
infinitely divisible distributions and about their classical and free cumulants. 
In Section 2 we explain the choice of the model (i.e., of the distributions 
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P^ and L^). In Section 3 we construct the distributions P^. Finally, the 
convergence in probability of the spectral distribution of a random matrix 
with distribution P^ to ^{fi) is proved in two steps. In the first one, we show 
the convergence when the Levy measure has compact support, and in the 
second one (in Section 6), we extend this result using approximation and 
compound Poisson distributions. The first step is achieved with the moment 
method, and is divided into two steps: convergence of the mean of every 
moment in Section 4, almost sure convergence in Section 5. The distributions 
L^ are constructed in Section 7.1, the convergence in probability of the 
symmetrization of the spectral distribution of |Mrf|, when M^ is distributed 
according to L^, is also divided in two steps. 

1. Preliminary results about infinitely divisible distributions. 



1.1. Definitions and the bijection ^. The results of this section concern- 
ing classical probabilities are in Gnedenko and Kolmogorov (1954) and in 
Petrov (1995); the results concerning free probabilities are in Bercovici and 
Voiculescu (1993) and in Bercovici, Pata and Biane (1999), except the con- 
tinuity of the inverse of the bijection ^, which is shown in Barndorff-Nielsen 
and Thorbj0rnsen (2002). A probability measure // on M is said to be *- 
infinitely divisible (resp. S-infinitely divisible) if for every n > 1, there exists 
a probability measure /u„ on M such that /u*" (resp. /U^") equals fi, which 
is equivalent to the existence of a sequence {fin) of probability measures, of 
a sequence {kn) of integers tending to infinity, such that /i*'^" (resp. fi^'^") 
tends weakly to fi. 

We can characterize *-infinitely divisible distributions (resp. ffl-infinitely 
divisible distributions) with their Fourier transform (resp. their Voiculescu 
transform). A probability measure ^ on M is *-infinitely divisible (resp. ffl- 
infinitely divisible) if and only if there exists a real 7 and a positive finite 
measure G on M such that its Fourier transform fl (resp. its Voiculescu 
transform ip^) has the form 



m 



exp< i-^t + 



MGI 



Atu 



1 



itu 
l + n2 



dG{u] 



(1) 



resp. ip^{: 



^7 + 



-^ for u= 

1 + tz 
^ z-t 



dG{t) 



In this case, the pair (7,G) is unique, and we denote fi = vl' (resp. v^ ). 
Remark. There exists other parametrizations of *-infinitely divisible 



G({0}),L(^)=/^i±^ 



dG(n) 



distributions: for example, denoting 7' = 7, a 

for all Borel set A of M \ {0}, one has /ugiR\|o|(l -^ u^)AL(u) < co, and 

m = exp(i7't - ^ + L6M\|0}(e^*" - 1 - w) d^(^))- 
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We now give the definition of ^, referred to in the Introduction. 

Theorem 1.1 (Bercovici-Pata's bijection). We endow the set of positive 
finite measures on M with the weak topology; the subsets {^-infinitely divisi- 
ble distributions} and {^-infinitely divisible distributions} are also endowed 
with the weak topology. 

1. The maps 

M X {positive finite measures} -^ {^-infinitely divisible distributions}, 

and 
M X {positive finite measures} -^ {^-infinitely divisible distributions}, 

are homeomorphisms and we have 

2. Let us define the map ^ , from the set of * -infinitely divisible distributions 
to the set of '^-infinitely divisible distributions, which maps, for all (7, G), 
the measure v2' to the measure v'^ . Then 

(a) ^ is an homeomorphism called Bercovici-Pata' s bijection, 

(b) for all fijV *-infinitely divisible distributions, ^{fi*^) = \I'(^)ffl^(i^), 

(c) Dirac measures are invariant under ^ : ^(^a) = ^a, 

(d) '^{N{m, r^)) is the semi-circle distribution with mean m and variance 
r^, which is Wm,2r{x)dx, with 

1 
2^r2 

(e) ^ , restricted to the Cauchy type, is the identity: for alia > 0, ^(Ca) = 
Ca, where Ca = ^J^, 

(f) for all sequence {^n) of probability measures on M, for all sequence 
{kn) of integers tending to infinity, the sequence ^"^^ converges weakly 
to a * -infinitely divisible distribution fi if and only if fi^ '^ converges 
weakly to ^(/i). 

Remark 1.2. In the texL the positive finite measure G is called the 
Levy measure of z/J' and v^ . We will use the two following properties: 

1. If the Levy measure of a ffl-infinitely divisible distribution u has compact 
support, then so does v [see Hiai and Petz (2000)]. 

2. i/g is symmetric if and only if v^' is symmetric, if and only if G is 
symmetric and 7 = 0. 



Wm,2r{x) = ^r^i^r'^ -{x- mff^ l\^_^\<2r, 
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1.2. Classical compound Poisson distributions, approximation of * -infinitely 
divisible distributions by ^-infinitely divisible distributions with compactly 
supported Levy measures. 

Definition 1.3. Let A be a nonnegative real, /? be a probability measure 
on M. Then the sequence of probability measures on M 

A\ A \ *" 

I--U0 + -P , n>l, 

n) n / 

converges weakly to a distribution noted vr* ^ , with Fourier transform 

75,(t) = exp(A(p(t)-l)), 
where p is the Fourier transform of p. 

Remark 1.4. tt*^ is z^?*^, with 

y_2 f U 

G = X— — i^dp{u), 7 = A/ — — i^<^p{u). 

We introduce now the compactly supported approximations of the positive 
finite measure G. 

Definition 1.5. Let, for G positive finite measure on M, t > 0, G^,Gt 
be the positive finite measures on M defined by 

G^,{A) = GiAn[-t,t]), GtiA) = G{A\[-t,t]) 

for all Borel set A of M. 

We define At > 0, the probability measure pt on M, and at £M with 

Xt= 5— dG(n), Pt = -, 5— dGt(u), 

JueR\[-t,t] u^ At u^ 



at = - {l/u)dG{u). 

Ju£M.\[-t,t] 



We will use the following approximation: 



/o\ w-t^ n 7,0 7+at.G'? 

(2) Vt>0 !//' = z^* *7r 



because one observes that '^*p x ~ ^"' with 
2 



H = \t-— — 5- dpf (u) = Gt, a = Xt -— — T,dpt{u) = -at. 

1 + n^ JneR 1 + n^ 
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1.3. Partitions, moments and cumulants of infinitely divisible distribu- 
tions. For every probability measure //, we will denote, when it is defined, 
by rnn{n) the nth moment of fi, which is J x^dfi{x). In this case, we will 
denote by <tn{fJ-) [resp. ^{fi)] its nth classical (resp. free) cumulant. Recall 
that [see Section 4 of Speicher (1994) or Section 2.5 of Hiai and Petz (2000)] 

(3) rnk{fi)= Y. n '^\V\ilJ') > 

7rePart(fc) V&tt 

denoted by CttCm) 

(4) mfc(^)= Y. I[^\V\it^) > 

7rGNC(fc) Vew 

denoted by ^^{fj.) 

where Part(/c) denotes the set of the partitions of {1,...,A;} and NC(A;) 
denotes the set of noncrossing partitions oi [k] = {1, . . . ,k} (a noncrossing 
partition of a finite totally ordered set / is a partition vr of / such that there 
does not exist x<y<z<t€l with x and z belonging to the same class 
and y and t belonging to another class). 

We will need the following proposition [part of which was proved in 
BarndorfF-Nielsen and Thorbj0rnsen (2004), but the proof we give here is 
shorter] : 

Theorem 1.6. Let n be a *-infinitely divisible distribution with com- 
pactly supported Levy measure, and let, for n integer, fin be a probability 
measure such that /i*" = /U. Then for each k > 1, the sequence {n x mk{fin))n 
tends to (J^fc(/u), which is equal to .^^.(^(yu)). 

Proof. By (3), one has 
nxmk{lJLn)=n ^ n^|V^|(^")= Yl ?^^~''''e:,r(/i) = (2^fc(/i) + o(l). 

7rePart(fc)ye7r'' ' 7rePart(fe) 

e;iv|{/i)/n 

Let us denote Vn =/^n"'- By part 2.(f) of Theorem 1.1, the sequence {vn) 
converges weakly to ^(/^). By Holder and Minkowski inequalities in tracial 
noncommutative VF* -probability spaces, every moment of v^ is bounded 
uniformly in n, so the cumulants of Vn tend to the cumulants of ^(;u). But 

by (4), 

7reNC{fc) Vevr' ' 7rGNC{fc) 



which tends to 

7reNC(fc) 



n 
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2. Free convolution and random matrices, choice of the models. For v 

probability measure on M, denote by P, the symmetrization of v^ which is 
the probability measure defined by viS) = \{y{B) + ^{—B)) for all Borel 
set B. 

For M Hermitian matrix, we will denote by fiM its spectral distribution, 
that is, the uniform distribution on its spectrum (with multiplicity). 

For M complex (possibly non-Hermitian) matrix, denote by j1\m\ the sym- 
metrization of the spectral measure of \M\, where |M| = \^M*M is the 
unique Hermitian nonnegative matrix such that M can be written M = 
U\M\, with U unitary. 

If M is a random matrix, fiM is a random probability mesure on the 
real line. For {Md)d sequence of random matrices, we will use the notion 
of convergence in probability for the sequence (/Ua/^) of random probability 
measures. 

The rest of this section may be skipped by the reader who wants to go 
straight to the result. We will only explain the choice of the models, that is, 
is of the family's P^ and L^ of distributions. 

Let us now explain in detail the choice of the family of the distributions 
P^, the distributions of the random Hermitian matrices. We would not go 
into as much detail for the distributions L^, which we construct in a similar 
way. 

The following theorem is proved in Voiculescu (1991) and in Pastur and 
Vasilchuk (2000) under more restrictive hypothesis, which can easily be re- 
moved using functional calculus. 

Theorem 2.1. Let n he a positive integer. Let ^i, . . . ,/i„, he prohability 
measures on M. Let, for d£ N*, (M^ )j=i^,,,^„ be a family of independent 
d X d Hermitian random matrices. We suppose that for all i = 1,. . . ,n, the 

(i) 

distribution of M^ is invariant under the unitary group^ s action, and fJ-^Ai) 

d 

converges in probability, when d^oo, to /ij. Then the spectral distribution 
of Y^^=i ^d converges in probability, when d^ oo, to ^ui ffl • • • ffl /U„. 

Let us consider a sequence (fin) of probability measures on M and a se- 
quence (kn) of integers tending to +00 such that /x*^" converges weakly to 

a probability measure fi on ]R. Let, for n G N, d G N*, (Af^*^)i<j<fc^ be a fam- 
ily of independent copies of a random Hermitian d x d matrix M^^m whose 
distribution is unitarily invariant. For every n G N, we suppose that l^M^^ 
converges in probability, when d — > 00, to /i„. 

Then we know that, for every n G N, the spectral distribution of Yl,i=i ^dli 
converges in probability, when d — > 00, to /^„ . 

Let us suppose that, on the other hand, for every d G N*, '}Zi=i^dn 
converges in distribution, when n ^ c«, to a random matrix M^. 
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We know, by Theorem 1.1, that /z^'^" converges, when n ^ oo, to the 
image ^(//) of /u by Bercovici-Pata's bijection. 

A natural question is the following: does the spectral distribution of M^ 
converge in probability, when d^ oo, to ^I'(/z)? 

In other words, is the limit, when d^ oo, of the spectral distribution of 
the limit, when n — > oo, of J2i=i ^dn ^Q^^^l to the limit, when n — > oo, of 

the limit, when d — > oo, of the spectral distribution of Yl,i=i ^dii^ 

The answer of this question is affirmative in our model [M^^n = U x 
diag(X„_i, . . . ,Xn^d)U* , U unitary Haar-distributed, independent of the i.i.d. 
random variables Xn^i, . . . , Xn^d with distribution fin]- It can be summarized 
in the following diagram: 

^il! + • • • + Afit"^ "^°" ' Pw 

d goes to OO d goes to oo 

i i 

spectral law: n^oo spectral law: 

The choice of this model is supported by the three following remarks: 

1. For d, n > 1, if /i„ = - ., /!^_f ^ , the expectation of the spectral distribu- 
tionofEr=iM«isi^. 

2. For any fixed d > 1, X)i=i-^rfn converges in distribution, when n ^ oo, 
to a distribution P^ which depends only on /u = lim„^oo /^n'^" • 

3. For every pair (//, v) of *-infinitely divisible distributions, similarly to the 
relation 

we have, for every d > 1, 

This property (and the fact that all formulas depend analytically on d, 
so could be extended to noninteger d) opens the perspective of a contin- 
uum between the classical convolution * and the free convolution ffl for 
infinitely divisible measures. 

Let us now explain how to construct the distributions L^. The following 
theorem is easily obtained combining the results of Haagerup and Larsen 
(2000) and Hiai and Petz (2000), and using functional calculus. 

Theorem 2.2. Let n he a positive integer. Let fii, . . . ,fin be probability 
measures on M. Let, for d> 1, (-/Vf^ )j=i^...^„ be a family of random d x d 
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(i) 

matrices with every M^ having a distribution invariant under the left and 
right actions of the unitary group. We suppose that, for every i = 1, . . . ,n, 

(i) 

the distribution of M^ is invariant under the left and right unitary group's 
actions, and that the symmetrization fi Ji) of the spectral distribution of 

\M^ I converges in probability to /ij. 

Then the symmetrization of the spectral distribution of 

n 

converges in probability, when d tends to infinity, to /Ui ffl • • • ffl /U„. 

Let us then consider, for // symmetric *-infinitely divisible distribution, 
a sequence (^„) of symmetric distributions and a sequence (A;„) of integers 
which tends to infinity such that /i„ " converges weakly to /U. Let d be a 
positive integer. If for all n, (M^*^)j=i^...^fc^ is a family of independent copies 
of C/diag(Xi, . . . , Xd)V , where U, V, Xi, . . . , X^ are independent, U and V 
are unitary Haar-distributed, and Xi,...,Xd are distributed according to 
/i„, then it appears that 

converges in distribution to a distribution L^ which depends only on /i. 

We will show that if M^ is distributed according to L^', then /i|Md| con- 
verges in probability to ^(/u). 

3. The distributions P^. E denotes expectation. For any distribution P 
and any function / on a set of matrices, Ep(/(M)) denotes //(M)dP(M). 
Tr denotes the trace. 

Theorem 3.1. Let ^i be an ^-infinitely divisible distribution. Let (fin) be 
a sequence of probability measures onM and (kn) a sequence of integers tend- 
ing to +00 such that the sequence /i*'^" converges weakly to /i. Let, for d > 1 
and n > 1, Q^" be the distribution o/C/diag(X„^i, . . . ,Xn^d)U* , where U is a 
Haar-distributed unitary random matrix, independent of the fin- distributed 
i.i.d. random variables Xn^i, ■ ■ ■ ,Xn^d- 

Then the sequence ((Q^")*'^") of probability measures on the space ofdx d 
Hermitian matrices converges weakly to a distribution P^. 

Moreover, Fourier transform of the distribution P^ on the space of d x d 
Hermitian matrices with the scalar product (Af, A^) 1— > Tr MN is given by this 
formula: for every Hermitian matrix A, 

(5) Ep^(exp(iTr^M)) = exp(E(dx V;,(('u,^u)))), 
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where 

• tp^ is the Levy exponent of fi, that is, the unique continuous function f 
from R into C such that /(O) = and the Fourier transform of fi is exp of , 

• (•,•) is the usual Hermitian product ofC^, 

• u = (ui, . . . , Ud) is a uniformly distributed random vector on the unit sphere 
ofC". 



It appears clearly that, for /x, u *-infinitely divisible distributions, P|^ 

d d 



Proof. We will show the pointwise convergence of the Fourier trans- 
form of the distribution (Q|^")*'^'" on the space oi dx d Hermitian matrices. 
Let Ahe a. dx d Hermitian matrix with spectrum a^W^. Let F„ (resp. F) 
be the Fourier transform of /x® (resp. /x® ). Then, when n tends to infinity, 
kn(Fn — 1) converges (uniformly on every compact set of M ) to the Levy 
exponent ijj oi ^®'^ (i-e., to V'u '^^ where 'ip^ is the Levy exponent of fx). 

We have 

E(«T,Mn-).fc„(exp(iTr AM)) = (EQMn(exp(iTr AM)))''". 



^d 



Recall Q^" is invariant under the unitary action, so 

E(Q^")*fen(exp(iTr AM)) = (E^Mn (exp(iTr(diag(a)M))))^". 

Q^" is the distribution of [/diag(X„^i, . . . ,X„ (i)C/*, where C/ is a Haar dis- 
tributed unitary matrix, independent of the /i„-distributed i.i.d. random 
variables Xn^i, • • • , ^n,d- So 

E(Qf"i)*fen (exp(i Tr AM)) 

= (E(exp(iTr(diag(a)[/diag(X„,i, . . . ,X„,d)C/* ))))'=" 

= (Ef expfz ^ akXn,i\uk,i\'^ \ j j 



i^n 



E LF„ Y,ak\uk,i\ 



\ \ \\k=i / ie[d]/ 

which can be written 



l + ^ElknlFnllj2^k\uk,in 



(recah [d] = {1, . . . ,d}). 
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But kn{Fn — 1) converges uniformly on every compact set to ip when 
n ^ oo , so we have 

E(QM.).fe4exp(iTr^M))"=^exp(EUM ^afc|nfc,/n J J J. 



It implies that (Q^")* " converges in distribution to a probability measure 
^ and that the Fourier transform of P^, c 
matrix A with spectrum a G W^, is given by 



P^ and that the Fourier transform of P^, evaluated on a d x d Hermitian 



E^.iexp{iTiAM))=expUUllj2''k\uk,in Jj )• 

But 'ip = 'il)®''-, so 

(6) EpM(exp(iTr^M))=exp(E(dx VM((^,a)))), 

d 

where (•, •) is the usual scalar product of W^ and Z = (|tiip, . . . , l^idP)) with 
u = (til, . . . , Ud) a uniformly distributed random vector on the unit sphere of 

Recall that the distribution of u is invariant under the unitary action, so 

E(d X V^((Z,a))) =E((i X ^^{{u,Au))). U 

Remark 3.2 (The Poisson case). One can already identify P^ when 
II = V{X) is the classical Poisson distribution with parameter A (denoted 
7r| ^ in Section 1.2). It is easy, using Fourier transform, to see that, in this 
case, P[^ is the distribution of 

X{dX) 
Yl Ud{k)ud{ky, 
k=l 

where {U(i{k))k>i is an independent family of uniformly distributed random 
vectors on the unit sphere of C^, independent of the P(dA)-random vari- 
able X{dX). 

Explicit computation of the Fourier transform o/P^ — the Gaussian case. 
In this section we give the distribution, the moments and the Fourier trans- 
form of the random variable Z appearing in (6) of the Fourier transform of 
P[^. In the following, we will only need the moments of Z. 

Proposition 3.3. Letu be a random vector of the unit sphere ofC with 
uniform distribution. Then the distribution of Z = {\ui\'^,. . . ,\ud\'^) on the 
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standard d-symplexe is the uniform distribution, that is, for every bounded 
Borel function f , 

E(/(|ni|2,...,|n,|2)) 

= {d-iy. / ••• / ' / j;i,...,Xd-i,l-Vxi dx. 

Jxi=0Jx2=0 Jxd-i=0 \ j^^ / 

To prove it, write u as a renornialized Gaussian standard vector on C^, 
and do an appropriate change of variables. 
We deduce, by induction on d, the fohowing: 

Proposition 3.4. For d>l, for a £ N*^, denoting s = !]«««) 

Eflm |2"i • • • luJ^"'') = (d- IV n»=i(«»0 < / n. (^-1)! 
H\ui\ \ud\ ) [ct ^J-(^^^_i), ^i^-J (^^^_i)r 

Remark 3.5. When fi = N{0, 1), that is, ■0/^(0 = ~V' Proposition 3.4 
ahows us to compute the Fourier transform. It appears then that, when 
Md has distribution P^, Md has the distribution of Nd + ,| X/^, where 

Nd G GUE{d,-g—j) [GUE{d,a'^) is the Euchdean space of Hermitian d x d 
matrices endowed with the standard Gaussian distribution with variance a'^] 
and X is a real standard Gaussian random variable independent of Nd- 

Proposition 3.3 allows us also to compute, by induction on d, the Fourier 
transform of the random variable Z. 

Proposition 3.6. Let d>2 be an integer and let a£W^ be such that 
the ak are pairwise distinct. Then 



ii^y^^yyi\u,^/jj — \ 


j=i llfc=i, 


..j,...,d^(«fc-"i)' 


This proposition, together with the formula 




EpM(exp(iTr^M)) 

d 




= exp{i7Tr(A)+d / 


E(e*"<^'">)-1- 


znTr(^) ■ 

d{l + u^)\ 


n2 



dG(u) 



E(<Z,a)2) 
2 U=0 

gives us the explicit computation of the Fourier transform of Pj. 
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4. Convergence of the fcth moment of the mean spectral distribution to 
the fcth moment of '^'(/x) when the Levy measure has compact support. 

4.1. Statement of the result, preliminaries for the proof 

Proposition 4.1. Let n be an *-infinitely divisible distribution with 
compactly supported Levy measure {in the sense of the definition given at 
Remark 1.2). Then we have 

Notation and preliminaries. Let, for n G N*, /x„ be a probability measure 
on M such that /_f*" = /i. Let us consider, for d > 1 and n > 1, (-^^„)i<j<n 
i.i.d. random matrices with distribution Q^". Then we know by Theorem 3.1 
that, for (i > 1, the sum of the M]*^'s (i = 1, . . . , n) converges in distribution 
to P(^ when n ^ oo. We know, by Theorem 1.6, that, for all /c G N*, the 
sequence nmk{l^n) is bounded, and so 

Vfc>l,Vci>l IEp^QTrM^)=JimE(iTr((|:M«) )). 

Let us then fix k £W. 

4.2. Computation of ¥,f,t^{2Ti M^) and proof of Proposition A.l. Let us 
define, for d, n > 1, 



a.,„ = E(±TV((5:M(J 



, i=l 



We have 



«^,n=XE( E n<i''^^' 




fe r=l 



/eN 

We will transform this sum by summing on the partitions. 

We denote by Bij(/) the set of permutations of a set /. Consider a par- 
tition TT of [n] (we have defined [n] = {1, . . . ,n}) and A: G [n]. We denote by 
7r(A;) the index of the class of k, after having ordered the classes according to 
the order of their first element [e.g., 7r(l) = 1; 7r(2) = 1 if 1 ~ 2 and 7r(2) = 2 
if 1 o*^ 2]. We denote, for l,n nonnegative integers, by Aj^, the number of 
one-to-one maps from [/] to [n], that is, n{n — 1) • • • (n — / + 1). 

The following lemma will be used quite often in the text. 
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Lemma 4.2. Consider k,n G N*. Consider (p: [n]^ -^ C such that 
V/G[n]^V(TGBij([n]) </,(a o /)= </,(/). 

Then 

/e[n]'= 7rGPart(fc) 

By this lemma, we have 

k 
i _ _/ ^r^ ,1^1 

ad. 

n \ 

TePart(fc) 

k 



\7rePart(fc) r=l / 

\7reNC(fc) r=l / " \7rePart(fc) r=l / 



denoted by ?;d,n 



denoted by w^ , 



Lemma 4.3. iei ir be a partition of a totally ordered finite set I. Then 
the following assertions are equivalent: 

(i) vr is noncrossing, 

(ii) there exists a class V of tt which is an interval, and vr \ {V} is a 
non- crossing partition of I\V . 

Using several times Lemma 4.3 and integrating successively with respect 
to the different independent random variables, we have 

^d,n = \Tv( E A\:\\{m\y\{^Jin)■IA 
V7reNC(fe) VGtt / 

7reNC(fc) ^->^-^ V^evr 

n — voo 

By Theorem 1.6, for every /c > 1, one has 

lim n X rukifin) =^k{^{fJ-))- 

n— »oo 

So for every d, 

}i^Vd,n= E n^|V'|(*(^))="^fc(*(/«))- 

7reNC(fc) FStt 
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To treat the term Wd,n, let us expand the trace: 

7r6Part{A,-) je[d]'= ^''=1 / 

7r^NC{fc) jk+i-=ji 

7r6Part{fc) rePart(fc) \r=l / 

7r^NC{fc) 

where for each r G Part(A;), r(/c + 1) = t(1). 

Using the fact that (M^*^)i<j<„ are independent copies of a matrix with 
distribution Q(^", we deduce 

^d,n = ^ E 4:' E 4' n %" f n ^^r(rM.+i)) 

7rePart(A:) rePart(fe) VGtt \reV J 

7r^NC(fc) 

=^ E 4^' E ^'n^f E nKM^^.(.+i),/.^n,u), 

7rePart(A:) rePart(fc) VStt \iG[d]^rey / 

7r^NC(fe) 

where [/ E U^ is Haar-distributed and independent of (X„,i, . . . ,^n,d)- So, 
applying Lemma 4.2, 

-.,n = ^ E 4-' E 4' 
7r6Part{A;) TePart(fc) 

7r^NC(fc) 

ySTT \crePart(y) rSV / 

integrating with respect to the X„^;'s, 

-.,. = ^ E 4^' E 4'n E 4^^' 

7rePart(fc) rePart(fc) FStt CTePart(y) 

7r^NC(fc) 

X Ef Jl ('Ur(r),a(r)^*r(r+l),<7{r)) I H "^bl (A*^) 



denoted by ad.T-,< 



2^ nl-l ^ "^ 

7rePart(A:) ^.^-^ rePart(fc) 
7r^NC(fc) "^5°1 



1 y^ An ■ Y^ ^|t| 
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vevro-ePartfy) n^ooJCTi ""ecr '^ « ' 

by Theorem 1.6 

when n — > cxo, for every vr ^ NC(A;), for every V £ tt, the only remaining 
aGPart(y) is {!/}. 
So one has 

7rePart(fc)TGPart(fc) V&tt \reV / 

7r^NC(fc) 

where u = {ui, . . . , u^) is a uniformly distributed random vector of the unit 
sphere of C^. 

Using the invariance of the distribution of u under the action of diagonal 
unitary matrices, one sees that for all k,l>0, z G [d] ,j G [d] , if 



\r=l r=l / 



then k = I and there exists a permutation (/) of [k] such that for all r, i^ 
So the preceding formula can be written 



d \d 

= ^ E E 4UMl[^\y\{f,)E(Y[u,^r)Urir+iX 

7r6Part(fc) Teacc(7r) VGtt \reV / 

7r^NC(fc) 

where for any finite totally ordered set / (in which the following element of 
any element x < max I is denoted by x + 1 and max/ + 1 = mini), for any 
partition vr of I, acc(7r) is defined to be the set of ir-acceptable partitions, 
which is the set of partitions r of / such that 

Vy G 7r,3 (/. G Bij(F), Vr G F T{r) = T{cl){r) + 1). 

Lemma 4.4. Let I be a finite totally ordered set, tt,t be partitions of I 
such that: 

• vr has a crossing (i.e., vr is not noncrossing) , 

• T is TT -acceptable. 
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Then we have 

(7) K| + |r| <|/|. 

Proof. We prove the lemma by induction on the cardinality of / (which 
is not less than four because tt has a crossing). 

• If the cardinality of / is four, then we can suppose / = [4] . We have 
vr = {{1, 3}, {2, 4}} and the inequality (7) is easy to verify because there are 
only three vr-acceptable partitions: 

{[4]}, {{1,2},{3,4}}, {{1,4},{2,3}}. 

• Suppose the inequality (7) proved when the cardinality of / is p, and 
consider / with cardinality p + 1, and vr, r partitions of / such that vr has a 
crossing and r is vr-acceptable. 

• If vr and r have no singleton class, then their cardinalities are not 
greater than |/|/2 and (7) is verified. 

• If vr has a singleton class {a}, then T{a) = T(a + l). This implies that 
if one removes the element a in /, the class {a} in vr and the element a of 
its class in r, then r stays vr-acceptable (and, clearly, vr keeps a crossing). 
So, by induction hypothesis, we have (|vr| — 1) + \t\ < \I\ — 1. 

• If r has a singleton class {6}, denote by V the class of 6 in vr and by 
(p the permutation of V such that for all r &V, T{r) = T{(j){r) + 1). We must 
have (j){b) -|- 1 = 6, so 6 — 1 ~ 6. Remove the element h in /, the class {6} in 
r and the element h oiV . Then, clearly, vr keeps a crossing. Define (j) to be 
the permutation of the "new" V by 

,/.(r), if0(r)/6, 

5-1, if0(r)=6. 

Then for all r in the "new" V , r and (j){r) are in the same class of the "new" r. 
It implies that r stays vr-acceptable. So, by the induction hypothesis, we have 

|vr|-F(|T| -1)<|/| -1. D 

Now recall Proposition 3.4: for a G N*^, denoting s = Yli oti, 

E(|ni|2"i • • • |urf|2"<') = {d- 1)! nti(«»0 < (^,). id - 1)! 



(s + d-1)! -' ' {s + d-l)\ 

But for vr, r E Part(A;), with r vr-acceptable, for all V ^ tt, there exists a G N'^ 
such that J2i Oii = \V\ and 



IE n ^^W"T(r-+i) ) =IE(|iiiP"^ •••|iid| 

\r£V J 



2arf 



18 F. BENAYCH-GEORGES 

So, by Proposition 3.4 we have 



m jd-iy- 



7r6Part(fc) T6acc(7r) V^Gtt ^i ' ' 



7r^NC(fc) 

Let C be real such that 



Vd>l,VsG[fc] (,!rJl^^<C7d- 



We then have 



<^ E E rfl^ldl'^l|^.(A^)|CWd"'. 

7r6Part(fc) TGacc(7r) 
7r^NC{fc) 

But according to (7), for all vr ^ NC(A;) and for all r € acc(7r), we have 
|t| + |vr| — A; < 0, so Proposition 4.1 is shown. 

5. Convergence in probability of the spectral distribution to '^(/x) when 
the Levy measure has compact support. 

5.1. Statement of the result and preliminaries to the proof. 

Proposition 5.1. Let n be an ^-infinitely divisible distribution with 
compactly supported Levy measure {in the sense of the definition given at 
Remark 1.2). Then the spectral distribution of a random matrix with distri- 
bution P^ converges in probability to ^(/i) as d tends to infinity. 

Notation and preliminaries. We keep the notation and the objects in- 
troduced in Section 4.1. We consider a sequence (M^) of random matrices 
defined on the same probability space such that for all d, M^ has distribution 
P^, and we will prove the almost sure weak convergence of the spectral dis- 
tribution of Mrf to *(//). It implies Proposition 5.1. Since ^(/u) is determined 
by its moments, the weak convergence of any sequence of distributions to 
^(/i) is implied by the convergence of all moments to those of ^(/i). 
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Let us fix k>l. We will show that almost surely, 

Var denotes the variance. 

Recall that by Borel-Cantelli's lemma, a sequence (Id)dgN of square- 
integrable real random variables converges almost surely to a real I if X)d(^(^) " 
/)^ and Z]dVar(y^) are finite. 

But we know that EpM(i TrM'') -mfc(^(/i)) = O(^). So it suffices to show 
that 



J2 VarpM (-TvmA < cx) 



a 

We will show that VarpM(iTrM'^) = O(^) using the formula 

/I . \ / I / / _J^ , .^ \ 

(8) 



Varp^QlVM^-) = Ji^n^Var (i Tr (( g M«) )) 



denoted by V^.n 

5.2. Computation o/ VarpM(^Tr Af'^') and proof of Proposition 5.1. We 
have 

\/e[.„]2fe" \r=l /" \r=fc+l // 



Let us apply Lemma 4.2: 



K,.= E 4-iEffiTvnMi;«))fiTr n Mi:«) 

7rGPart{2fc) \\ r=l / \" r=fc+l 

E A'^^' Alr^lE f i Tr n Af£«)) IE f i Tr n m£ 

,ePart(fc) V" r=l / V" r=l 



^n ^n ^ ^ ^ -^^ 11 ^"-'d,n j ^ ^ ^ -^^ 11 ^■'■'. 

''ri,vr2SPart(fc) 

We split the sum into two parts: in the first one we sum over the par- 
titions of [2fc] which can be split into two partitions tti and 7r2 respectivly 
of {1, . . . , /c} and of {/c + 1, . . . , 2/c}, in the second one we sum over other 
partitions of [2A:], 

''ri,'''"26Part(fc) 
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xE 



i^E< 



(Mr)) 



r=l 



E 



i^IlK 



(772 (r)) 



r=l 



.r=l 



7rePart(2fc) 

Let us expand the trace: 

T/ , — _ \^ MI'^iI + I'^^I _ /l|vri| /l|7r2| 



1 



2k 



^Tr n<i^^MX n <i^^^ 



, r=fc+l 



d^ 



''ri,7r2€:Part(fc) 



4 E n(<^^")>,>..)^( E n(<n^^'")>,>...^ 

/ V ie[d]'= '"=1 > 



V je[d]fe r-=l 

jk+i-=ji 



3k+i--=n 



-i E 4-' E iE[(M£^^)),,,---(<'^"k.. 



d2 



7rePart(2A:) i6[d]2'= 



x(M^ 



(7r(fc+l)), 



.(7r(2fc)), 



'd,n >jk+i,ik+2 ^ d,n 'J2k,jk+i'' 

We apply Lemma 4.2 once more: 



d2 



'!'ii''r2GPart(A;) 
A; 



ri,r2£Part(fc) 



xE 



ll(^<i,n )ri(r),ri(r-+l) 



r=l 



E 



n«n 



(7r2(r))^ 



Lr=l 



/r2(r),T2(r+l) 



2fc 



ll(^rf,n. )r(r),T(r+l) 



r-=l 



+ ^ E 4^' E 4^'^ 

7rePart(2fc) rGPart{2fc) 

where for any partition r of [/c], r(A; + 1) denotes t(1), and for any partition 
r of [2k], l<r<2k + l,we define 

T(r), ifs^{A: + l,2fc + l}, 

f{r) = { t(1), if r = fe + l, 

t(/c+1), ifr = 2A: + l. 



r»- 



tiM" 



Since (M^ „)i<i<n ai'e independent copies of a matrix with distribution v^^ , 
we have 



(i2 



''ri!'r2SPart(fc) 



Ti,T2GPart(A:) 
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niE, 



VeTTi 



n ^ri(r),ri(r+l) 
r&V 



W^ 



J V£lT2 



+ i E 4^^' E 4'nE«3^^ 



7rGPart(2fc) TGPart(2A:) VStt 



n ^r2(r),r2(r+l) 
n ^^r(r),T(r+l) 

rev 



1 ^ (4^^1+1-1 -Ai-14'^^i) Y. 4-Ui-i 



(i2 



i"l,i"2GPart(fc) 



■'"i,T2SPart(fc) 






E n "T2(r-),«,^^r2(r+l),«,-'^n,i, 



E 11 '^T{7-),lrUf{r + l),lr^ri 

laUWrmv 



+ ^ E 4' E 4'niE 

7rePart{2fc) rGPart{2fc) VStt 

where C/ is a unitary Haar-distributed random matrix, independent of {X^^i^ ■ ■ ■ , Xn^d) 
So, after application of Lemma 4.2, 



d- 



T'i,'n'2€Pa.rt{k) 



Ti,r2£Part(fc) 



^X„ 



^ n ^( E 4 n '"n{r),(T(r)^ri(r+l),<T(r)^n,fT{r) 

VGtti \(TGPart(y) rSF / 

^ n^ E 4 n '"^2(r),<7(r)^r2(r+l),<7(r)^n,<7(r) 



ye7r2 \crGPart(y) r£V 



+ 



\t\ 



l E 4' E 4 

7r6Part(2fc) r6Part(2fc) 

ni^f E 4'n-^M,^ 

VStt \crePart(y) rev 



VStt \crePart(y) 

integrating with respect to the X„ ;'s. 



VnA 



1 

d2 



E 

''ri,7r2GPart(fc) 



4ki| + |7r2| <|7ri| ,|7r2| 



(r)tt^(j.+l)^o-(j.)X„^o-(.r) I, 



Inl Jral 



n 



|7ri| + |7r2| 



E 4"'^, 

Ti,T2ePart(fc) 
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X n E n^-I^I^^E 

FG7r2(TGPart{y) 



11 '^Ti{r),a-{r)'^Ti(r+l),a-{r) 
IreV 

n ■"r2{r),<7{r)^r2(r+l),<T{r) 

rev 



Y[nx m|^|(/i„,) 

uGcr 

J|n X m|j,|(/i„) 



. vGa 



+ i E 4 E 4'n E -^-'^'4'^' 



7rGPart(2fc) rGPart(2A;) VGtt o-ePart(y) 



3i<k<j,ir^j 
X E Jl 'Ur(r),(T(r)%(r+l),<T(r) 

-rev 
Let n tend to infinity: 



Y[nxm\^\{fin)- 



■uGo- 



1 



>4,d — ,2 Z^ 



(i2 



.kil + lTTzl _ .kil .Kal , , , , 



''ri,7i"2SPart(fe)'' 



n 



171-11 + 17121 



'ri,T2ePart(fc) 







xR E ^4^' 

- — ^ d' 



xE 



r&V 



i(r),(T(r)^ri(r+l),(T(r) 



Jl nx m|„|(/i„) 



. I'Gcr 



by Theorem 1.6 



n E iiiA 



i-kl a\^\ 



V&1T2 crePart(y) n— Kx> |o-| 
> 0^ 



xE 



11 ''^T2(r),(T(r)'^r2(r+l),cr(r) 

rey 



Jl n X m|„|(/i„) 



. "UGcr 



by Theorem 1.6 



+i E ^ E 4'n E ^4^^' 



7rePart(2fc) ^-«^— ' TePart(2fc) ye7rCTePart(y)n^oo |<t 

TT n — >CXZ) 1 

3i<k<j^i'^j — * 1 



xE 



11 ^T(r),(T(r)^f(r+l),o-(r) 
reF 



H^ nxm|^|(/x„) 



v£a 



by Theorem 1.6 
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when n tends to infinity, for every partition vr (or vri or 112), for every V £ tt, 
the only resting a E Part(y) is {V}. So one has 

Varp^f^TrM^' 

= i E E 4"' n ^^^1 VI (/^)IE ( n "r(0%(r+l) ) , 

7rGPart(2fc) rGPart{2fc) VStt \reV / 

where u = (ui, . . . , Ud) is a uniformly distributed random vector of the unit 
sphere of C . 

But recall that by invariance of the distribution of u under the action of 
diagonal unitary matrices, for all k,l>0, i£ [d]'',j G [d]', if 



\r=l r=\ ) 



then k = I and there exists a permutation </> of [/c] such that for all r, iy. 
So the preceding formula can be written 

VarpMl^TrM* 
<i\d 

7rePart(2fc) r6adm(7r) VStt VreV / 

where adm(7r) is defined in the following way (splitting the set [2A:] in two 
disjoint sets [/c], [2/c] \ [A;]): for any pair (/, J) of disjoint finite totally ordered 
sets, for any partition vr of / U J, adm(7r) is defined to be the set of vr- 
admissible partitions, which is the set of partitions r of / U J such that 

VT/GTT, 3(/)EBij(y),VrGF T(r) = r(0(r) + 1), 

where for any x € I (resp. x G J), x + 1 denotes the element following x in 
/ (resp. J). 

Lemma 5.2. Let (I, J) be a pair of disjoint finite totally ordered sets, 
TT,T partitions of I U J such that: 

• there exists i £ I,j £ J , with i ~ j , 

• T is TT -admissible. 

Then we have 

(9) |7r| + |T|<|I| + |J|. 
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This inequality can be proved by induction, the proof is analaguous to 
the one of (7). 

Recall (Proposition 3.4) that for a G N"^, using the notation s = J^i'^t^ 

But for every 7r,r G Part(2A;), with r vr-adniissible, for every V £ n, there 
exists a G N*^ such that J2i Oii = \V\ and 



Vrey / 



2a^l . . . Ui ,|2adi 
I "'rfl J 



So, by (10), we have 
VarpMl -TrAf 



d \d 



d^ A,/'' ll^|v|(/^)(H/|-J (|^M^_i)r 



7rePart(2A:) VGtt 

3i<k<j,ir^j 
r£adin(7r) 

Consider C < oo such that 



Vd>l,VsG[2/tl (g!) % ^^ /^\, <Cd- 
~ ^ ^ ^ ' {s + d- 1)1 ~ 



We have 



2k 



\yd\<^ E d\^\s^\\c^{f,)\c\-\d- 

7r6Part(2fc) 
TGadm(7r) 

But according to (9), for every vr G Part(2A;) such that there exists i <k < j 
with i ~ j and for every vr-adniissible r G Part(2fc), we have |r| + |7r| — 2A; < 0, 
so 

Va,>,.(ilVM'=)=0(-l 

and Proposition 5.1 is proved. 

5.3. Applications to GUE and sums of independent projections. This 
section is not necessary for the rest of the text. 

Proposition 5.1 contains the almost sure convergence of the spectral distri- 
bution of the matrices of GUE{d, ^xj-) ^° ^^^ semi-circle distribution, where 
GUE[d,a'^) is the Euclidean space of d x d Hermitian matrices with the 
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scalar product Tr(- x •), endowed with the standard Gaussian distribution 
with variance o"^. 

Indeed, let {Nd)deN* be a sequence of random matrices such that for 
every d, the distribution of A^^ is the one of a matrix of the GUE{d, -^j^) 
[we do not do any hypothesis about the joint distribution of (A^d)deN*]- Let 
X be a real Gaussian standard random variable, independent of (A^d)deN* ■ 
We have seen to Remark 3.5 that for d eW, M^ := N^ + -7= • Id has 

distribution P^ ' . We have proved that fiM^ converges almost surely to 
the centered semicircle distribution with variance 1. So /Utv^, which is equal to 
6_ X * fj,Md 1 converges almost surely to the centered semicircle distribution 

with variance 1. 

Another consequence of Proposition 5.1 is the following one. Recall that 
for all A > 0, the Marchenko-Pastur distribution with index A is the image, 
by the Bercovici-Pata bijection, of the classical Poisson distribution V{\) 
with index A. 

Proposition 5.3. Let, for all d>l, {ud{k))k>i he an independent fam- 
ily of uniformly distributed random vectors on the unit sphere of C^ . Then 
for all A > 0, the spectral distribution of 

d' 

^Ud{k)ud{k)* 

k=l 
converges in probability to the Marchenko-Pastur distribution with index A 
when d,d' tend to infinity and the ratio d'/d tends to A. 

The proof of this result, which uses tools introduced in the following 
section, is in the Appendix. 

6. Convergence in probability of the spectral distribution Md to *(/i) 
without condition on the Levy measure. 

6.1. Convergence in probability of a sequence of random distributions to 
a deterministic distribution. We will denote, for z G C, by SRz and Qz its 
real and imaginary parts. Let us define, for i^ probability measure on R, 

/^:C+ = {zGC;9z>0}^C, 

duiu) 



/neR U — Z 

Then f^ is a holomorphic function on C"*^, |/i/(-z)| < c^, and the map 
{probability measures on M} -^ M"*", 

(/xi, ^2 )^sup{ 1/^,(2;) - f^^2{z)\;Qz> 1} 
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is a distance which defines the weak topology. 

So, for (Mrf)rf>i sequence of Herniitian random matrices and p probabihty 
measure on M, we have equivalence between: 

(i) the spectral distribution of pM^ converges in probability to p, 
(ii) for every e > 0, 



P( sup 



^-l^{^,{Md))-fp{z] 



>e] — ^ 0, 



where, for M Hermitian matrix and z G C \ M, 9^z(M) = (M — z)~^ . 

6.2. Statement of the theorem and scheme of the proof. 

Theorem 6.1. Let p be an ^-infinitely divisible distribution. Let, for 
d > 1 , Md he a random matrix with distribution P^ . 

Then the spectral distribution pM^ of M^ converges in probability to ^(p). 

Scheme of the proof: 

1. Notation, approximation of M^ by M^. 

2. Upper bound, for a > 0, of p(rg(A'^^) > da) uniformly in d > 1. 

3. Conclusion. 

6.3. Proof of Theorem 6.1. 

6.3.1. Notation, approximation of Md by M^. Let 7, G be such that p = h'2' • 
Recall [equation (2)] that for t > 0, denoting: 

1. By Gf and Gt the positive finite measures on M, 

G^,{A) = G{An[-t,t]), Gt{A)=G{A\[-t,t]) 

for all Borel set ^ of M. 

2. The at the number -L6M\[_t,f] :^dG(u). 

6. By pt,nut the measures i^* ,1/* , 

we have the following: 

(i) p = pt^vt, so for every d, Ma has the distribution of M^ + Nj^, where 
MJ and Nj^ are independent random matrices with respective distributions 
P^* andP^S 

(ii) Vf is the weak limit, when n — > co, of 

1 ]So-\ Pt 

n J n 



FREE INFINITELY DIVISIBLE DISTRIBUTIONS 27 

with 

1 + ii^ .^. X 1 l + n2 



dG(n), Pt = -; 5— dGj(u). 

/«GK\[-t,t] "" Ai U^ 

So for all d > 1 , the distribution P^* of A'^^ is the weak limit of the distribution 

of Sr=i -^dn ' where, for every n > 1, {N^ J^ )i<i<n are independent copies 
of C/diag(X„,i,...,X„,d);7* with: 

(a) {Xn^i, ■ ■ ■ ,Xn^d) i-i-d. random variables with distribution 

1 K>oH Pt, 

n ) n 

(b) \J unitary Haar-distributed random matrix, independent of (X„^i , . . . , X^^d) 

6.3.2. Upper bound, for a > 0, ofP{rg{N^) > da) uniformly in d>l. We 
denote by rg(M) the rank of a matrix Af . 
Let a be a positive real. 
Rank is a lower semi-continuous function, so 



E(rg(iV*))<Jiin^Efrgf5:iV*;^'^ 

<„iii^]E('x:rg(ivi?) 



n 



lim^E(rg(iV*;W)) 



n— >oo 

= hm nE(rgC/diag(X„,i,...,X„,d)f/*) 

d 

= hm n5]p(X„,,^0) 

n— >cxD -^ — ' ' 

1=1 

= lim ndp(X„,i/0) 
= lim nd — - 

n — *oo 11 

= dXt. 
So we have 

E(rg(iV*))<dAi. 
We deduce, with the Chebyshev inequality, that, for every a > 0, 

(11) P(rg(iVi) > da) < -LE(rg(iV*)) < ^. 

da a 
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6.3.3. Conclusion. Let e, r] be positive reals. Let us show that there exists 
an integer do such that, for every integer d>do, 



Pi sup 



ilV(lH,(Mrf))-^(^)(z) 



> e] <!]. 



Choice oft>0. When t tends to +00, the real at tends to and the pos- 
itive finite measure G^ converges weakly to G. So, by Theorem 1.1, z^^ *' * 
converges weakly to u^ . In other words, ^(/xt) converges weakly to ^(fi). 
So there exists Ti > such that, for all t>Ti, 

(12) sup I /vi,(^,,) {z) - U(^^) {z)\<-. 

Oz>l "J 

When t tends to +00, the real At tends to 0, so there exists T2 > such 
that, for every t > T2, 

er] 



(13) 

Let t = max(Ti , T2) . 

For every d> 1 , we have 



Ai< 



12 



Pi sup 



-T¥(^,(Md))-^(^)(z] 



> e 



Pi sup 



;Tr(lK,(Mi + iVi))-/*(^)(z) 



>e 



and 



(14) 



TTi9\,{M', + N',))-hUz) 



< 



d 
+ 



Tr(lK,(M* + iV*)-lK,(M*)) 



-Tr(lK,(M*))-/*(^,)(. 



d 



+ \h(^lt)iz)- f^!(^,){z)\. 



Let us deal with the first term of the sum (14): 

We know that for every complex dx d matrix M, |^ TrM| < — g- 



^rg(M), 
where ||M|| is the operator norm of M associated to the canonical Hermitian 



normonC^, and ||^^(M* + iV*)-9^,(M*)|| < ||^^(M* + iV*)|| + ||9^,(M*)|| < 

isz — 

Moreover, for all pair Ad,N of Hermitian matrices, for all z G C \ M, 
9i,{M + N)-m,{M) = -9\,{M + N)N9\,{M). So rg(9^,(M* + iV*)-9^,(M*)) < 
rg(iV*). 

So 

(15) 



;Tr(lH,(M* + iV*)-lH,(M*)) 



2 

d 



<3i'g(^i)' 
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but for all d > 1 
1 



Pi sup 



d 



TT{^,{Mi + Ni)-^,{Mi)) 



> 



<p 



rg(ivi: 



< 



< 



p 

e 
V 



d 

rg(iVd 



> 



*^,>^^ 



6 J 



by (11) 
by (13). 



By inequality (12), the third term of the sum (14) is < | as soon as Qz > 1. 

Let us now deal with the second term of the sum (14). The Levy measure 
of fit (in the sense of the definition given at Remark 1.2), which is G^, is 
compactly supported. By Proposition 5.1 and by the other results of Section 
6.1, there exists an integer dg such that, for every d> do, 



Pi sup 



-Tr(lH,(Mi))-/*(^,)(. 



d 



- 3J 2 



Then for all d> do, replacing the terms of the sum (14) by the upper 
bounds we just gave, we have 



Pi sup 



TV(^,(Mrf))-^(^)(z 



>.l<| + | + 0. 



So, we have 



lim P sup 



Tr(lH,(Mrf))-/*(^)(^ 



>e 



and Theorem 6.1 is proved. 

7. Study of the non-Hermitian modeL 

7.1. The distributions L^. This section is the analogue, for non-Hermitian 
matrices, of Section 3. The distributions L(^ are defined by the following the- 
orem, the proof of which is analoguous to the one of Theorem 3.1 using the 
polar decomposition of non-Hermitian matrices and the bi-unitarily invari- 
ance of the distributions K.^" . 

Theorem 7.1. Let fi be an ^-infinitely divisible distribution. Let (//„) be 
a sequence of probability measures on M and (kn) be a sequence of integers 
tending to -|-oo such that the sequence /i"^" converges weakly to fi. Let, for 
d>l and n > 1, K^" be the distribution o/C/Diag(X„^i, . . . ,Xn^d)y, where 
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U, V are independent unitary Haar- distributed dx d random matrices, inde- 
pendent of the fj.n- distributed i.i.d. random variables X„,,i, . . . ^X^^d- 

Then the sequence {{K^")*^") of probability measures on the space ofdx d 
complex matrices converges weakly to a distribution L^. 

Moreover, the Fourier transform ofL,'^ on the Euclidean space of complex 
d X d matrices endowed with the scalar product {M,N) i— > ?R.{Tr M* N) is 
given by the following formula: for all complex dx d matrix A, 

(16) Elm (exp(iK(TrA*X)))=exp(E(dx^^(SR((n,Af ))))), 

where : 

• ip^ is the Levy exponent of /i, 

• (•,•) is the canonical Hermitian product of C , 

• u = {ui, . . . ,U(i), V = {vi, . . . ,Vd) are independent random vectors, uni- 
formly distributed on the unit sphere of C . 

Remark 7.2. 1. Notice L^ * L;^ = h'^*" . 

2. When ^ = A^(0, 1), L^ is the distribution of a matrix [Mij] with (3f?Mjj,9Mjj)i<ij<(i 
A^(0, ^)-distributed i.i.d. random variables. 

3. The same construction can be done with rectangular bi-unitarily invari- 
ant random matrices. It leads, when the dimensions of the matrices tend 
to infinity in a certain ratio, to probability measures which are infinitely 
divisible with respect to a certain convolution. The studying of this convo- 
lution has led the author to construct a new noncommutative probability 
theory, called the rectangular free probability theory, which allows us 
to understand the asymptotic behavior of rectangular random matrices, 
as free probability theory describes the asymptotic behavior of square 
random matrices. It might give rise to a publication. 

7.2. Convergence of the kth moment to the kth moment of ^(/i) when 
the Levy measure is compactly supported. The purpose of this section is to 
show the following result: 

Proposition 7.3. Let n be a symmetric *-infinitely divisible distribu- 
tion with compactly supported Levy measure. Then for all integer k, 

lEL^("^fc(A|A.f|)) - "2fc(*(Ai)) = 0(^- 

ProOF. First, for every complex d x d matrix M, for all integer k, 
mk{iJ'\M\) is null if k is odd and is equal to ^Tr(MM*) '^ if k is even. 
As jJL is symmetric, ^(/i) is symmetric. So it suffices to show that, for all 
keW, 



ELMQTr(MAr)^) -m2fc(^(^)) = oQ). 
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Let, for n G N*, /in be the probability measure such that //*" = /i. Consider, 
for d> 1 and n > 1, (-^^„)i<i<n i-i-d. random matrices with distribution 
K^". Then we know by Theorem 7.1, that, for every d>l, the sum of the 
M^*^'s (z = 1, . . . ,n) converges in distribution to L^ when n goes to oo. 

We know, by Theorem 1.6, that, for ah k £N* , the sequence (n x mk{Hn))n 
is bounded, and so that, for ah fc,(i G N*, 

(17) IE(m2.(/i|M|))=Jim IE(iTr(^(^X:^/«) (E^^?)) )• 

Let us fix k gW . We are going to use (17). 
Let, for d, n > 1, 

/l //" \/" vvi-v 



6.,„ = lE(^^Tr((^MW)(^M^: 
We have 






^^,n=^Tr(E( j: n 



, ,{/{2r-l)) .M{2r))* 



VG{l,...,n}2 

Let us transform this sum using partitions (Lemma 4.2). Moreover, from 

(i) 

now on, we do not write anymore the index d in M^^, 



bd,n = ^ Tr ( E ( Y. 4"Im('^(i))m^'^(2))*^(.(3)) . . . m(-(2'=))* I ) . 

V \7rGPart{2fc) // 

But 

E(m(i)*m(i) • • • m(i)*) = E(M^i)M^i)* . . . m(i)) = = ma+i(/i„)/d, 

^ V ' " V ' 

21+1 alterned factors 21+1 altcrned factors 

E(Af^^)*M(i) • • • M«) = E(M^i)M^i)* • • • Mi^>) = m^iifin)^- 

21 altcrned factors 21 alterned factors 

So, for vr G NC(2A;), using many times Lemma 4.3 and integrating succes- 
sively with respect to the different independent random variables, we obtain 

^TV(E(mWi))*mW2))mW3))* ...j\^W2fc)))) = ^^(^^)_ 
Proceeding then like in Section 4.2, we show easily that 

IEL^("^2fc(A|Af|)) -"l2fc(^(/i)) 

= 1 E ^l['^'''''^'r(/^) n ^( n Ur{r)Vr{r+l) Yl ^T(r+l)'^T{r) ) , 

7r,TGPart(2fc) VGtt \ reV reV ) 

7r^NC(2fc) ^ °dd 
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with r(2A; + 1) = r(l) and where u = (ui, . . . ,Ud),v = (fi, . . . ^Vd) are inde- 
pendent uniformly distributed random vectors of the unit sphere of C^. But 
as we have ah'eady seen, by invariance of the distribution of u under the ac- 
tion of unitary diagonal matrices, for every pair (7r,r) of partitions of [2A;], 
if 

n^( n Wr(r)Wr(r+l) H ^r(r+l)^r(r) J 
Vevr V rev r&V / 

r odd r even 

is nonzero, then for every class V of vr, there exists (j), permutation of V , 
which maps odd numbers to even numbers and vice versa, such that for all 
r S y, T{r) = T{(p{r) + 1). It implies that r is vr-acceptable. Using inequality 
\t\ + |vr| < 2k [equation (7)], the inequality on the moments of a uniform 
random vector on the sphere of C^ (Proposition 3.4), we deduce, as in Section 
4.2, that 



7r,TePart(2fc) V€w \ r<=V reV 



{r+l)VT{r) ] =0[ - 



T,TePart(2fc) 
7r^NC(2fc) 



r odd 



So 



VA;GN E^.{m,{fl\M\))-m,{^{fi)) = 0[^,. □ 



7.3. Convergence in probability to ^(/u) when the Levy measure has com- 
pact support. The purpose of this section is to show the following result: 

Proposition 7.4. Let ^ be a symmetric ^-infinitely divisible distribu- 
tion with compactly supported Levy measure {in the sense of the definition 
given at Remark 1.2). Let, for each d, M^ be a random matrix with distri- 
bution L^. 

Then the symmetrization of the spectral distribution of |Mrf| converges 
weakly to ^(/i) when d goes to infinity. 

Proof. We will show inequalities that would imply almost sure con- 
vergence of the symmetrization of the spectral distribution of \Md\ to ^(^) 
if the matrices M^ {d>l) were defined on the same probability space. So 
let us suppose that the matrices are defined on the same probability space. 
We keep the notation and objects introduced in Section 7.2. Since *(^) is 
symmetric and determined by its moments, the convergence of a sequence 
of symmetric distributions to ^(/i) is implied by the convergence of all the 
moments of even order to those of ^(/u). 
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Let us fix k>l. We will show that almost surely, 



^TriM*,M,)^''^m2ki^ifi)). 



But we know that 



ELMQTr(MM*)^) -m2fc(^(^)) = oQ 



So it suffices to show that 



Var,.QTY(MM*)'=)=o(l 



We will do it using the formula 

Vari^.(iTr(MiVr)'=) = Jl^ Var (^ixr (^(^X:<i) (E^^J 
Proceeding like in Section 5.2, we obtain 
VarQTV(M^*Mrf)^- 

7r,T6Part(4fc) Ven \ reV reV J 

-, .^„, . .IT . r odd r even 

where u = {ui, . . . , U(i),v = {vi, . . . , v^) are uniformly distributed independent 
random vectors of the unit sphere of C^. 

But for all couple (7r,r) of partitions of [4A;], if 

n^( n '"r(r)'^f(r+l) H ^r(r)%(r+l) 1 
r odd ^ even 

is nonzero, then for all class V of vr, there exists a permutation <^ of V, 
which maps even numbers to odd numbers and vice versa, such that for all 
r €V, T{r) = f{(p{r) + 1), which implies that r is vr-admissible. Using the 
inequality |r| + |7r| < 4A: [equation (9)] and the inequality on the moments of 
a uniform random vector on the sphere of C'' (Proposition 3.4), we deduce, 
as in Section 5, that VarLM(^Tr(MM*)'=) = O(^). D 

Remark 7.5. In the case where /z = A^(0, 1), we have a new proof of a 
well-known result: the spectral distribution of a Wishart dx d matrix with d 
degrees of freedom converges almost surely, when d tends to infinity, to the 
distribution of X^ when X is a centered semi-circular random variable with 
variance 1, which is [see Speicher (1999)] the Marchenko-Pastur distribution 
with parameter 1. 
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7.4. Convergence in probability of ji\mA ^^ ^(/i) in the general case. 

Theorem 7.6. Let ji be a symmetric *-infinitely divisible distribution. 
Let, for d>\, M^ be a random matrix with distribution L^. 

Then the symmetrization ft\MJ of the spectral distribution of |M(^| con- 
verges in probability to ^(^). 

The proof is quite similar to the one of Theorem 6.1. 

Proof of Theorem 7.6. 

Notation, approximation of M^ by M\. Let G be the symmetric positive 
finite measure on M such that /i = z^*' . Recall [equation (2)] that, for t > 0, 
if: 

1. Gj and Gt are the positive finite measures 

GO(^) = G(^n[-i,t]), Gt{A)=G{A\[-t,t]) 



for all Borel set A of 

O.G? 
lJ't = l^* 

then we have 



o o.G? , o,Gt 

2. ^t = i^* and ft = I/*' , 



(i) ji = fit * i^t, so for all d> 1, Ma has the same distribution as M^ + 
A'^^, where Mj and A'^^ are independent random matrices with respective 



distributions L^* et L^*, 



(ii) vt is the weak limit, when n — > oo, of 

1 OoH Pt 

n J n 



for 



■dG(n), Pt = \ 5— dGj(u). 



MeK\[-t,t] u^ Ai u^ 

So for all (i > 1, the distribution L^' of N\ is the weak limit of the distribution 
of Z]r=i -^d'n ' where for all n > 1, {N^^)i<i<n are independent copies of 
^diag(X„,i, . . .,XnA)V with: 






(a) (Xn,i, . . . , Xn^) i.i.d. random variables with distribution (1 — -^)6q + 

(b) [/, T^ unitary Haar-distributed random matrices, independent of {Xn^\, . . . , Xn^d) 
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In the same way as in Section 6.3.2, we show that, for all a > 0, 



(18) 



P(rg(iV*) > da) < 



Xt 



Let us denote, for p probability measure on M, p^ the distribution of X'^ 
when X is a random variable with distribution p. 

Consider e,r] > 0. 

Let us show that there exists an integer do such that, for all d>do, 



P( sup 

.Sz>l 



TV(9^,(M,*M,))-/^(^)2(z 



> e <r]. 



Choice oft>0. When t tends to +oo, the measure G^ converges weakly 

G^* n C 

to G. So, by Theorem 1.1, u^ ' converges weakly to u^ . In other words, 

*(/xt) converges weakly to ^(/-f). So '^{pt)^ converges weakly to '^{pf'. 

Hence, there exists Ti > such that, for all t>Ti, 



(19) 



sup \hipt)2{z) - ^(^)2(Z)| < -. 



When t tends to +oo, the real At tends to 0, so there exists T2 > such 
that, for all t > T2, 



(20) 



^*-24- 



Let t = max(Ti ,T2). 
For all d > 1 , we have 

1 



P( sup 



-T\-(lK,(M,*M,))-^(^p(z) 



> s 



P( sup 



Tr(9^,((Mi* + N*f)iM', + iVi))) - /*(^)(z) 



>e 



Hence, 



- Tr(^,((Mr + iVr )(M^ + K))) - /*(^)(z) 



< 



- Tr(9^,((M** + iV**)(M* + iV*)) - 1H,(M**M*)) 



d 



TV(fn,(MfMi))-/^(^,)2(z) 



+ \h{^lt)^{z)-U{^^r{z)\■ 



But for all pair {M,N) of Hermitian matrices, for all z G C \ M, 9^2 (M 
N) - ^z{M) = -^,{M + N)N9\,{M). 
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(21) 



So denoting A*^ = MI^*NI^ + N^/{M^a + iV^), we have 

1 
d 
< 



Tr(lH,((AC + N^niM'd + N'd))) " h^i^ 



- TV(^,((M** + iV**)(Mi + iVi))A^9^,(M**M*)) 



+ \h{titriz)- f^,^^)2{z)\. 



The conchision is similar to the one of Section 6.3.3. D 

APPENDIX 

We prove Proposition 5.3. Denote J2k=i '^d{k)ud{k)* = N^^d' ■ As explained 
in Section 6.1, it is equivalent to prove that, for every e > 0, 



P( sup 

,S2>1 



TV(5H,(iVd,dO)-/. 



*{^(A)) 



(^) 



d,d'—too „ 

>e I >0. 

d'/d~\ 



By Proposition 5.1, ^('P(A)) is the limit of the spectral distribution of a 
random matrix with d 
is the distribution of 



random matrix with distribution F^ ■ But, as noticed in Remark 3.2, P^ 



X(dA) 

Md:= J2 Mk)ud{kT, 
fc=i 

where X{dX) is a 'P((iA)-random variable, independent of the sequence {ud{k)). 
So it suffices to prove, that for every e > 0, 



P( sup 



Tii^^iN, 



dA' 



^y^ziMd)) 



d,d'— >oo „ 

>e\ >0. 

d'/d~\ 



But it was noticed in Section 6.3.3, equation (15), that 



\TT{^,{Nd,d:)-^,{Md)) 
a 



<-vg{Nd,d'-Md), 
a 



which is not greater than ||X((iA) — d'\, which converges in probability to 
zero, by the weak law of large numbers. 
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